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ABSTRACT 

The boundaries of the Uranian e, a, and (3 rings can be fitted by Keplerian 
ellipses. The pair of ellipses that outline a given ring share a common line 
of apsides. Apse alignment is surprising because the quadrupole moment of 
Uranus induces differential precession. We propose that rigid precession is 
maintained by a balance of forces due to ring self-gravity, planetary oblateness, 
and interparticle collisions. Collisional impulses play an especially dramatic role 
near ring edges. Pressure-induced accelerations are maximal near edges because 
there (1) velocity dispersions are enhanced by resonant satellite perturbations, 
and (2) the surface density declines steeply. Remarkably, collisional forces felt 
by material in the last ~100 m of a ~10 km wide ring can increase equilibrium 
masses up to a factor of ~100. New ring surface densities are derived which 
accord with Voyager radio measurements. In contrast to previous models, 
collisionally modified self-gravity appears to allow for both negative and 
positive eccentricity gradients; why all narrow planetary rings exhibit positive 
eccentricity gradients remains an open question. 



Subject headings: celestial mechanics — planets and satellites: individual 
(Uranus, rings e, a, and f3) 



1. INTRODUCTION 

Each narrow eccentric ring surrounding Uranus is composed of particles moving on 
nested elliptical orbits. The outer and inner edges of a given ring define ellipses having 
semi-major axes a ± Aa/2 and eccentricities e ± Ae/2, where Aa -C a, Ae <C e, and e < 1. 
Observed values of a, e, Aa, and Ae for the Uranian e, a, and (3 rings are listed in Table |I[ 

Remarkably, the set of ellipses describing an individual ring share a common line of 
apsides. Apse alignment is surprising because the oblateness of Uranus causes orbits of 
particles with different semi-major axes to precess differentially. Timescales for differential 
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precession in the absence of other forces are extremely short; in the case of the Uranian e 
ring, the inner edge would precess a full revolution relative to the outer edge in 175 years. 
Rigid precession of an eccentric planetary ring has remained a problem in ring dynamics for 
over 20 years. 

Goldreich & Tremaine (1979, hereafter GT) proposed that apse alignment is maintained 
by self-gravity. Their theory predicts that the eccentricity gradient across the ring, 

_ de 
qe = a d-a> 

must be positive. A positive eccentricity gradient in an apse-aligned ring implies that the 
ring is narrowest at periapse and widest at apoapse. Gravitational forces between particles 
are therefore greatest near periapse. Material in the inner half of the ring pulls radially 
inward on the outer half at periapse, generating a differential precession which exactly 
cancels that due to planetary oblateness. 

Though the prediction that q e > accords with observations of all known narrow 
eccentric rings, the standard self-gravity model (hereafter SSG) predicts Uranian ring 
masses that are too low compared to those inferred from Voyager radio occultations. Ring 
masses based on observations exceed predictions by factors of at least ~3 (e ring) to ~50 
(a and (3 rings) (Tyler et al. 1986; Gresh 1990; see also the reviews by Esposito et al. 1991 
and French et al. 1991). Low surface densities are particularly problematic for the a and 
(3 rings. With SSG surface densities, torques exerted by inner shepherd satellites would be 
insufficiently strong to confine the a and f3 rings against drag from the distended exosphere 
of Uranus (Goldreich & Porco 1987, hereafter GP). In addition, as discussed by Graps et al. 
(1995), shapes of the e ring surface density profiles as derived from occultation light curves 
do not accord with SSG predictions. 

This paper points the way towards resolving these problems. In §||, a theory of 
collisionally modified self-gravity (hereafter CMSG) is qualitatively described. A simple 
quantitative model is set forth in §||], in which new surface density profiles are derived for 
the e and a rings that are in better agreement with observations. In §[|, implications of 
our solutions for torque balance, the role of planetary oblateness, and the value of q e are 
discussed. Directions for future research are summarized in §^|. 
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2. QUALITATIVE SOLUTION 



For simplicity, consider an apse-aligned eccentric ring having constant, positive q e 
across its width. The ring is filled with spherical particles of internal mass density p and 
radius r, and the ring surface density is given by S. Let n and Co be the mean motion and 
apsidal angle, respectively, of a ring particle. Subscripts i and b denote quantities evaluated 
in the ring interior and near the ring boundary, respectively. Variables subscripted with 
p or s are associated with the central planet or shepherd satellite, respectively, and take 
their usual meanings. The dimensionless strength of the quadrupole moment of the planet 
is given by J 2 . Numerical estimates in this section are made using parameters appropriate 
for the e ring. 

A key ingredient missing in the SSG model is an accounting for interparticle collisions. 
Since ring optical depths r measured normal to the orbital plane are typically of order 
unity, each particle collides, on average, a few times with its neighbors every orbital period. 
Only modest collisional impulses per unit mass and time, of order 0.1 cms -1 per orbit, are 
required to generate differential precession rates comparable to those induced by planetary 
oblateness (GT). Velocity dispersions of order q ~ 0.1 cms -1 in the ring interior are not 
unreasonable: both the Keplerian shearing velocity across a particle diameter, ~ 3nr, and 
the escape velocity from the particle surface, ~ r^8Gp, are of that order for the meter-sized 
bodies that plausibly compose the ring. 

Although a single collision can impart an impulse of dynamically significant magnitude, 
multiple collisions experienced by a particle in the ring interior leave its precession rate 
largely unaltered. A particle in the ring interior is struck by its inside neighbors about 
as frequently and as forcefully as by its outside ones. Differential precession across 
the ring induced by smooth internal pressure gradients occurs on timescales of order 
27rEnae/| VP| ~ 27rnaeAa/cf ~ 10 6 (0.1 cms" 1 /^) 2 yr, much longer than misalignment 
timescales set by planetary oblateness (cf. GT). Here the height-integrated pressure 
P ~ Scf is taken to vary over a lengthscale Aa. 

Conditions are dramatically different near ring edges. Pressure-induced accelerations 
are maximal there because (1) velocity dispersions are enhanced by resonant satellite 
perturbations, and (2) the surface density declines steeply (Borderies, Goldreich, & 
Tremaine 1982). The velocity dispersion near the ring boundary could be as high as 




(2) 



where d ~ 10 3 km is the ring-satellite separation, and w r is the width of the annulus perturbed 
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by the satellite. To order-of-magnitude, the latter is given by w r ~ a ^JM S /M P ~ 1km, 
the distance from the resonant edge at which nested periodic orbits cross. Equation (0) 
is derived by equating the rate of energy dissipation by collisions in the perturbed zone, 
~ 7iT,aw r clnr, to the rate of energy deposition by the satellite, ~ 3nTd/2a ~ 9nT l nac*d/2, 
where T is the satellite-induced confining torque whose magnitude equals that of the viscous 
torque, ~ 37r£c 2 a 2 , in steady-state. 

A particle on the ring edge experiences a radially directed, collisional acceleration 



C ~ — — ~±-^-r~±Q,nf, (3) 

where the upper (lower) sign applies to the outer (inner) ring edge. Here P ~ S c 2 is taken 
to vary over a radial lengthscale, A, of order the local ring thickness, c^jn. In a q e > ring, 
collision rates are highest near periapse. At the periapsis of a ring boundary, the radial 
acceleration, C, generates a differential precession rate, A(du/dt)c ~ —C/nae, relative to 
the precession rate at the ring midline. This collision-induced rate is greater than the local 
differential rate due to planetary oblateness, A(du/dt) , by a substantial factor: 



A(du/dt) c TCb/ae ^ c h ^ 



A(doo/dt}o T(2lAa/8a)J 2 n( J R p /a) 2 lcms" 

Self-gravity maintains apse alignment against differential precession caused by 
planetary oblateness and interparticle collisions. For self-gravity to enforce rigid precession 
near ring edges, surface densities there must be higher than those predicted by SSG. At 
ring boundaries, self-gravitational attraction must balance the extra repulsive acceleration 
due to collisions. To estimate the surface density near the edge, £&, equate the collisional 
acceleration, C ~ c 2 /A, to the gravitational acceleration from a wire of linear mass density 
£&A located a distance A away: 



c 2 b /\ ~ 2GS 6 . (5) 

Take A = c^/n ~ 50 m to obtain 

E b ~c 6 n/2G~10 3 - Cb gem" 2 , (6) 
lcrns -1 



which is greater than corresponding SSG predictions by factors ^ 40. Equation (||) is 
equivalent to the condition that Toomre's Q be of order unity at the edge. 
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These endwires of mass ~ 2iraT lb \ = ncla/G constitute new boundary conditions not 
found in SSG. Gravitational forces from massive endwires induce substantial differential 
precession in the ring interior. For self-gravity to maintain apse alignment in the interior, 
surface densities there must also be greater than those predicted by SSG. 



Divide the region occupied by an apse-aligned, constant q e = aAe/Aa ring into an even 
number N of equally spaced intervals. The center of the j th interval contains an elliptical 
wire having mass rnj, semi-major axis aj = a + [j — (N + l)/2] Aa/N, and eccentricity 
Cj = e + [j — (N + l)/2] Ae/N. Denote by Aj(du/dt) the precession rate of the j th wire 
relative to the precession rate of a test particle at the ring midline. Uniform precession 



Subscripts O, G, and C denote contributions from planetary oblateness, self-gravity, and 
interparticle collisions, respectively. The first two terms are given by 



3. QUANTITATIVE MODEL 



requires 






(8) 



and 





where 




(cf. GT). 

For Aj{duj/dt)c, the following simplistic prescription is adopted: 
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,d(jQ s 



+ [q e H(q 2 e ) cl/Xnae] (1 - x/X) if x = (j - 1/2) Aa/N < X 



A,(— > c =<j -[q e H{q 2 e )c 2 b /Xnae)(l-y/X) ify = (N — j + 1/2) Aa/N < X (10) 
otherwise. 

Thus, Aj(du/dt)c is non-zero only within intervals 5a = X from each edge; there, its 
magnitude rises linearly from to q e H(q 2 )c 2 / Xnae. As a crude justification for this 
maximum value, approximate the collisional acceleration as 



C « -— 

w ±17^ n f ( n ) 

A(l - g e cos/) 

where the upper (lower) sign applies to the outer (inner) edge. Here the pressure gradient 
is taken to vary inversely as the separation between streamlines. Insert C into Gauss's 
perturbation equation for duj/dt and average over true anomaly: 



dt nnae Jo 

= TqeH(q 2 e )c 2 b /Xnae. (12) 

For the constant q e ring models presented here, 

X = c b /n, (13) 
so that the only remaining free parameter is c&. Equation ([H|) is relaxed for § [4.2.2| and 



4.3.2. 



Note that this prescription for Aj{du/dt)c ignores the decrease in velocity dispersion 
from Cb at the ring edge to q in the ring interior. The decline in velocity dispersion occurs 
over a length scale of order w r . This length scale is large compared to Q,/n so that the 
gradient of velocity dispersion does not give rise to a significant radial acceleration. 

For a given value of q,, equations (|7|), (]|), (§), (|I0"D , and fll3|) comprise N equations in 
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N unknowns Solutions for surface density profiles at quadrature for various values 

of Cft are displayed in Figure 1, for parameters appropriate to the e and a rings; models 
for the /3 ring are nearly identical to those of the a ring. In CMSG models, higher surface 
densities near ring edges are evident, as are higher total ring masses. 



4. DISCUSSION 

4.1. Surface Density Profiles and Torque Balance 

Simple CMSG models, while not fully realistic, demonstrate the existence of a new 
class of self-gravity solution, that obtained by accounting for the modification of ring 
boundary conditions by interparticle collisions. Remarkably, forces felt by material in the 
last ~100 m of a ~10 km wide ring can increase equilibrium masses by factors up to 100. 

Large S-band opacities measured by Voyager, which are incompatible with SSG surface 
densities (see, e.g., the review by Esposito et al. 1991), can be reconciled with average 
CMSG surface densities of ~75-100gcm -2 for the e, a, and (3 rings. Moreover, CMSG 
models predict that surface densities near ring edges are higher than those in the interior. 
This behavior is reminiscent of the "double-dip" structure seen in occultation light curves 
for the e and a rings (see, e.g., the review by French et al. 1991). 

Greater ring masses as implied by CMSG resolve problems associated with exospheric 
drag that were pointed out by GP for rings a and j3. For the remainder of this subsection, 
numerical estimates will be made for the a ring; similar conclusions hold for the j3 
and e rings. Surface densities are scaled to a typical CMSG value in the ring interior 
of £ = 75 gem -2 . An inner shepherd satellite exerts a repulsive, non-linear torque at 
first-order Lindblad resonances of magnitude 



L ^ l0p s R 3 s Z 2 n 2 a 7 17 / g V f R s \ 3 / p s \ / 500km \ 

nl ~ M%d l^gcm" 2 / UOknJ l^l^gcm" 3 /^ d ) QIg ' 

(14) 

where the satellite radius, R s , is scaled to the Voyager upper limit of 10 km (Smith et al. 
1986). The shepherding torque exceeds the magnitude of the drag torque exerted by the 
Uranian exosphere, 



1 Reflection symmetry about the ring midline reduces the number of equations necessary to N/2. Typically 
N > 2000 wires are needed to converge to within 10% of the solution for N — > oo. 
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Fig. 1. — Surface density profiles at quadrature for rings e and a in CMSG (q, ^ 0, solid line) 
vs. SSG (cf, = 0, dashed line) models. When the total CMSG ring mass, M CMS g, greatly 
exceeds the total SSG ring mass, Mgsc, it is found empirically that Mcmsg oc c^/A 1 / 2 oc cj 2 . 
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T d « -Anm H n H v T na 3 Aa « -4 x 10 16 f — " H „ j erg . (15) 

VIO 3 cm"v 

Here n# = 7 x 10~ 6 e 32 - 4R p/ a cm" 3 is the number density of hydrogen atoms of mass 
in the exosphere, and vt ~ lkms" 1 is their thermal speed normal to the orbital plane 
(Broadfoot et al. 1986). That > \T d \ ensures that the inner shepherd prevents ring 
particles from spiraling in towards Uranus. 

Estimates of viscous torques T v also require revision. For a ring undergoing Keplerian 
shear, with minimum kinematic viscosity n(E/p) 2 , the viscous torque is given by 

_ 37m 2 £ 3 a 2 18 / E \ 3 /l.5gcm- 3 \ 2 , , 

T v ~ 5 ^2.5xl0 18 - — 2 erg 16 

p l \75gcm l ) \ p J 

(GP). That T v 3> \T d \ ensures that ring particles on the outer edge press against the inner 
Lindblad resonance established by the outer shepherd. 

Conclusions drawn from comparisons between and T v are on less sure footing. For 
the choice of scaling parameters, the latter exceeds the former, contrary to the requirement 
of the standard theory of shepherding that the torques be equal. This might be construed 
as evidence that the angular momentum luminosity in the ring interior is reduced below T v 
by the non-Keplerian shear associated with a non-zero q e (Borderies, Goldreich, & Tremaine 
1982; GP). However, the numerical estimates for the two torques differ only by a factor of 
a few. The shepherding torque should be evaluated using surface densities near the edge, 
which CMSG predicts are higher than those in the interior; this would increase the estimate 
of T,j. Uncertainties in the choice of parameters preclude drawing any conclusion other 
than that these torques are of the same order of magnitude. 

4.2. Relative Importance of Planetary Oblateness 

4.2.1. J 2 = vs. J 2 + 

What does CMSG predict if J 2 = 0? Figure 2a displays the answer for the e ring, for 
0, = 2 and 3 cms" 1 . In contrast to SSG, a non- vanishing equilibrium surface density does 
not require a finite planetary oblateness; self-gravity can be balanced entirely by collisional 
pressure gradients. For the a and (3 rings, solutions with and without J 2 are practically 
indistinguishable for c& > 0.5 cms -1 . The influence of J 2 on the equilibrium solution 
diminishes as Aa decreases or as c& increases. 
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Fig. 2. — (a) CMSG e ring models for which J 2 is reduced from its nominal value of 3.35 x 1CT 3 
(Elliot & Nicholson 1984) (solid line) to (dotted line). As c& is increased from 2 cms -1 (lower 
two curves) to 3 cms -1 (upper two curves), the influence of J2 diminishes. For the a and j3 
rings, CMSG models with and without J2 are practically indistinguishable for c?, > 0.5 cms -1 
(data not shown), (b) CMSG e ring models for which q e = ± 0.626. Contrary to SSG models, 
a positive q e is not required to obtain an equilibrium solution. 
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4-2.2. Empirical Scaling Relations for J 2 = 

For J 2 = and fixed ring geometry, the surface density at quadrature near a given edge 
scales as 

Z b (0<\x\&X) = -^f(\x\/\), (17) 

where \x\ measures distance from the edge, q, and A are the same free parameters as in 
Equation (|10D, and / is a dimensionless function of the similarity variable |x|/A. Well away 
from ring edges, the surface density at quadrature scales as 

^(\x\>\) = QjL=g(\x\/Aa), (18) 
where g is another dimensionless function. The total ring mass scales as 



M~ ^^Aa/A. (19) 



4.3. Value of q e 

4.3.1. Sign of q e 

Figure 2b displays a CMSG model for the e ring obtained by reversing the sign of q e . 
In contrast to SSG, a positive eccentricity gradient is not necessary in CMSG to obtain 
an equilibrium solution. This resurrects the problem of why all known eccentric planetary 
rings, including the Titan and Huygens ringlets around Saturn, are narrowest at periapse 
and widest at apoapse. 

It is possible that equilibria obtained using q e < are unstable. To address this 
issue, a preliminary investigation of ring stability for an N = 4 ringlet model has been 
undertaken. Forces due to pressure gradients are included only for the first and fourth 
ringlets. Collisional accelerations are treated as if they arise from anti-self-gravity forces 
(self-gravity with the sign of the acceleration reversed); i.e., collisional shear stresses are 
ignored. In this crude approximation, equilibria are found to be stable regardless of the 
sign of q e ; small deviations from equilibrium masses result in apsidal librations (Borderies, 
Goldreich & Tremaine 1983). It remains to be seen whether collisional shear stresses alter 
stability properties. 
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Another possibility is that initial conditions set the sign of q e . If the ring were initially 
uniform in width as a function of azimuth, then planetary oblateness would determine the 
initial sense of differential precession within the ring. The resultant narrowing of the ring 
width near a true anomaly of / = — n/2 would cause a positive eccentricity gradient to 
grow by self-gravity. Under this hypothesis, an iV = 2 ringlet model incorporating forces 
from self-gravity and planetary oblateness yields the following time evolution for the apse 
and eccentricity differences between outer and inner ringlets: 



5uj = -A sin Oubt (20) 
5e = Ae(l - cosfi lib t) (21) 

where A > and Qm, are the amplitude and frequency, respectively, of libration (cf. 
Borderies, Goldreich, & Tremaine 1983). Note that the time-average of 5e is positive. 
Inelastic collisions would damp librations and the ring would eventually settle into an 
equilibrium for which q e > 0. 



4-3.2. Magnitude of q e Near Ring Boundaries 

It has been assumed that the eccentricity gradient, q e , is finite out to the last 
A = Cb/n ~ 50 meters of ring material. A finite q e is necessary to generate a non-zero 
azimuthal average of the collisional acceleration [see equation (|T2"D]. The simple quantitative 
model of employed the observed value of q e averaged over the entire ring width. The 
true value over the last few hundred meters of ring material is unknown. 

In the case of the best-studied e ring, Graps et al. (1995) combined Voyager 
photopolarimeter and radio occultation measurements to infer the eccentricity gradient as 
a function of semi-major axis. They found that q e decreases over the last ~5 km from its 
nearly constant value of ~0.65 in the interior to ~0.35 near the edge. The radial resolution 
of their study was between 1 and 2 km. 

A decrease in q e towards ring boundaries is theoretically plausible. Distortions in 
a circular ring can be described by the change in separation, 5r, between neighboring 
streamlines of the form 



5r oc cosm(0 — fi pat t) , (22) 
where f2 pat is the pattern speed of the distortion and m is an integer. A constant q e ring 
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that precesses rigidly in the quadrupole field of the central planet is equivalent to a distorted 
circular ring for which m = 1 and f2 pat = {dCb/dt)q. Resonant satellite perturbations, which 
enhance velocity dispersions within a distance w r ~ 1 km of ring edges, are characterized 
by much higher values of m = 2a /3d 3> 1 and f2 pat = Satellite-induced disturbances 
might therefore reduce the local value of q e . A decrease in q e over a distance w r near ring 
boundaries is roughly equivalent to setting A = w r in equation fllOD . By the scaling relations 
([UP and ([19|) , this would reduce surface densities and total ring masses shown in Figure [l] 



This work is primarily a demonstration that interparticle collisions near ring boundaries 
play a crucial role in determining ring masses under the self-gravity hypothesis. The nature 
of ring boundary conditions has not been calculated in rigorous detail; instead a prescription 
motivated by order-of-magnitude arguments is provided for collision-induced precession 
rates. Numerical simulations incorporating shepherd satellites will help to determine the 
actual 3-dimensional collisional stress tensor and eccentricity gradient everywhere within 
the ring. 

Why all narrow eccentric rings surrounding Uranus and Saturn are narrowest at 
periapse and widest at apoapse remains to be understood. Stability analyses incorporating 
collisional shear stresses may reveal that rings having q e < are unstable. Alternatively, 
the sign of q e may be set by initial conditions. Scenarios for ring formation — e.g., the 
catastrophic disruption of a small moon — require further elucidation. 

Viscous damping gives rise to small differences between apsidal angles of neighboring 
streamlines (Borderies, Goldreich, & Tremaine 1983). For a given apsidal shift of 5 Co <C 1, 
the difference between the azimuth of maximum streamline separation and the azimuth 
of apoapse is given by the "pinch angle", 5(f> = arctan(e&D/£e) ^> 5Co. The pinch angles 
calculated by Borderies et al. (1983) for their N = 2 streamline models of the Uranian 
and Saturnian ringlets are suspect, however, because they neglect the boundary effects 
highlighted in the present work. A careful calculation of 5<p{a) that incorporates viscous 
drag and the global effects of resonant forcing by shepherd satellites has yet to be performed. 
Upcoming observations of narrow Saturnian ringlets by the Cassini Orbiter might test the 
predictions of such a calculation, thereby furnishing a powerful diagnostic of stresses within 
ringlets. 



by a factor of Jw r n/cb ~ 4. 
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Table 1. Parameters of Eccentric Uranian Rings a 



Ring a(km) Aa(km) e(xl0 3 ) Ae(xl0 3 ) 



e 


51149 


58.1 


7.936 


0.711 


a 


44718 


7.15 


0.761 


0.076 


a 


45661 


8.15 


0.442 


0.066 



a Values taken from Tables I and VII of French et 
al. (1991). 



